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Abstract. We present results for the chemical evolution of the Milky Way including predictions for 
elements from De\ 
us to draw import; 
galaxy formation. 



elements from Deuterium to Europium. A comparison with the most accurate and recent data allows 
us to draw important conclusions on stellar nucleosynthesis processes as well as on mechanisms of 
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INTRODUCTION 

> , 

The Milky Way is the best studied system and a great deal of accurate data such as high 
resolution spectroscopy relative to the solar neighbourhood have recently accumulated. 
Detailed chemical evolution models able to follow in detail the evolution of the abun- 
dances of several chemical species in the interstellar gas are also available and by means 
of them we can compare theory with observations and infer many important constraints. 
In particular, by adopting different set of stellar yields for stars of all masses we can test 
the stellar nucleosynthesis and the origin of elements. Particularly important are the plots 
relative to abundance ratios, such as [X/Fe], as functions of stellar metallicity, measured 
by the quantity [Fe/H]. Recently, Cayrel et al. (2004) have provided a sample of abun- 
dances for extremely metal poor stars in the Galactic halo, in particular for stars with 
metallicity as low as [Fe/H]=-4.0dex. The elements studied are mainly a-elements (O, 
Mg, Si, Ca, S, Ti), Fe-peak elements and heavier.Abundances for C and N in extremely 
metal poor halo stars were also measured (Spite et al. 2005; Israelian et al. 2004). All of 
these data represent an invaluable template for testing chemical evolution models. Here 
we test the two-infall model (Chiappini et al. 1997) for the chemical evolution of the 
Milky Way on a large sample of data including the abundances discussed above plus 
the abundance trends recently derived for light elements such as deuterium. The main 
ingredients of the chemical evolution model are: i) the initial conditions, ii) the star for- 
mation history, iii) the initial mass function (IMF), iv) the stellar nucleosynthesis, v) the 
rate of gas accretion. In the two-infall model the main assumptions are that the halo and 
thick disk formed by gas accretion on a relatively short timescale (1-2 Gyr), whereas 
the thin disk formed by means of an independent accretion episode on longer timescales 
increasing with galactocentric distances. In particular, the timescale for the formation of 
the solar neighbourhood is 7-8 Gyr, as imposed by the stellar metallicity distribution. 
The assumed rate of gas accretion (primordial gas) is: 

IR=A(R)e- t/TH{R) +B(R)e^ t - t,m " ) ^ D ^ (1) 
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with the timescale for the first infall episode Th(R)= 2 Gyr for the halo, and the 
timescale for the disk infall Td(R) linearly increasing with R which implies an "Inside- 
Out Model". The time t max is the time of the maximum infall on the thin disk. 

The star formation rate is: 

SFR = vo&o$ u (2) 

with k\ = 1 .5, k2 = 0.5, and a tot being the total surface mass density and o gas the surface 
gas density. A threshold density for the SFR of lM & pc~ 2 is also assumed. This threshold 
regulates the star formation regime producing an oscillatory behaviour for the SFR when 
the gas density is close to its value (see Figure 1). 

The assumed IMF is multi-slope (xi, x%, ..) and follows the prescriptions of Scalo 
(1986). 



Element production 

The element production in stars can be summarized as follows: 

• Low & intermediate mass stars (0.8 < M/Mq < 8) produce mainly 4 He, 12 C, 14 /V 
and heavy s-process elements (Ba, Y, Sr) (in stars with masses from 1 to 3M & ). 

• Type II SNe (M > 8M©) produce mainly a-elements (O, Ne, Mg, Si, S, Ca), part 
of Fe and other Fe-peak elements as well as r-process elements (Eu, Ba). 

• Type la SNe produce mainly Fe-peak elements (~ 0.6 — 0JM & of 56 Fe per super- 
nova). 

• Novae can possibly be important in producing 7 Li, 13 C, l5 N, 11 0. 1 Li can be 
produced also in SNII, red giant, asymptotic giant branch stars and cosmic rays. 

• Deuterium is only destroyed in stars, whereas 3 He is destroyed but also produced 
in stars from 1 to 3M Q . 

Type la SN rates 

The type la SNe are the most important Fe producers in the universe and represent 
a fundamental ingredients in galactic chemical evolution models. The most popular 
scenarios for the progenitors of Type la SNe can be summarized as follows: 

• Single degenerate scenario. This is the classical scenario of Whelan and Iben 
(1973), namely C-deflagration in a C-0 white dwarf reaching the Chandrasekhar 
mass, Mch, after accreting material from a companion which can be either a red 
giant or a Main Sequence star.In this case the maximum mass allowed for both 
the primary and secondary star in the binary system is 8M and consequently the 
minimum timescale for the occurrence of type la SNe is tsNia min =0-03 Gyr (Greggio 
& Renzini 1983; Matteucci & Recchi, 2001). 

• Double Degenerate scenario. This is the merging of two C-0 white dwarfs, due 
to gravitational wave radiation, which explode by C-deflagration when the M C h is 
reached (Iben & Tutukov 1984)The minimum timescale for the occurrence of such 
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FIGURE 1. The star formation rate as predicted by the two-infall model by Chiappini et al. (1997) for 
the solar neighbourhood. Note the oscillatory behaviour at late times in the thin disk evolution which is a 
consequence of the assumed star formation threshold. 



systems is given by the lifetime of a 8M , as above, plus the gravitational time 
delay t S ma min = 0.03 + A% rav =0.03+ 0.15Gyr (see Tornambe & Matteucci, 1986). 
• A more recent model by Hachisu et al. (1999) is based on the single degenerate 
scenario but with a metallicity effect. In particular, no type la system can form 
if the metallicity of the progenitor stars is lower than [Fe/H] = —1.0. This model 
suggests also that only systems with secondary stars with masses lower than 2.6M 
can be accepted. This leads to a minimum time for the occurrence of type la SNe 
no shorter than tsma min — 0.33 Gyr, to which one should add the metallicity delay, 
namely the time taken by the chemical enrichment process to reach a metallicity of 
[Fe/H]=-1. 

It is important to know the minimum timescale for the occurrence of type la SNe as 
well as the time for the maximun in the type la SN rate, which depends not only upon 
the SN progenitor model but also upon the SFR. 



Results 

The G-dwarf metallicity distribution. In Figure 2 we present the predicted G-dwarf 
metallicity distribution compared with the observed distributions. The good fit to the 
data is obtained by means of a timescale for the formation of the thin disk in the solar 
vicinity of 7-8 Gyr. A much lower timescale would produce a very bad fit to the observed 
data since it would move the peak of the distribution toward too low metallicities (see 
Kotoneva et al. 2000). The stellar metallicity distribution is, in fact, very sensitive to the 
star formation history which in turn depends strongly on the rate of gas accretion. 
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FIGURE 2. The continuous line represents the predictions of the two-infall model with a timescale of 
disk formation of 8 Gyr.The histograms represent the observational data as indicated in the figure. 



Abundance ratios versus iron abundance. In Figures 3 and 4 we show the predic- 
tions relative to the [a/Fe] ratios versus [Fe/H] compared with a complete set of data 
including Cayrel et al.'s (2004) data, as described in Francois et al. (2004), who tested 
different sets of yields on the data. Their conclusion was that the yields of massive stars 
by Woosley & Weaver (1995, WW95) and those of Iwamoto et al. (1999, 199) for type 
la SNe are the best to reproduce the observations altough some variations are required 
in the yields of some elements such as Mg, Zn and some Fe-peak elements. The good 
fits shown in Figures 3 and 4 were obtained by including these variations. 

A remarkable feature of Figures 3 and 4 is the lack of spread at very low metallicity, 
suggesting that the halo of the Milky Way was already well mixed at [Fe/H]=-4.0. This 
should be taken into account in inhomogeneous models of chemical evolution of the 
halo. 

The evolution of light elements. In Figure 5 we present the evolution of deuterium 
during the lifetime of the Galaxy in the solar vicinity. The predictions are from the 
models of Chiappini et al. (2002) and Tosi (1988). The observed values at the time 
of birth of the solar system and in the local intertellar medium (ISM) are also shown. 
The assumed primordial value for D is the one suggested by the WMAP experiment 
(Spergel et al. 2003). The large spread shown for the local abundance of D, as measured 
by FUSE and HST-GHRS probably indicates the effects of D depletion into dust grains. 
It is important to notice that, in order to reproduce the data, the astration factor for D 
should have been lower than a factor of 1.5 (see Romano et al. 2003). 

Very heavy elements. The s- and r- process elements are generally produced by 
neutron capture on Fe seed nuclei. The former are formed during the He-burning phase 
both in low-intermediate and massive stars, whereas the latter occur in explosive events 
such as Type II SNe. Recently, Francois et al. (2006 and this conference) have measured 
the abundances of several very heavy elements (e.g. Ba and Eu) in extremely metal poor 
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FIGURE 3. Predicted and observed [a/Fe] ratios for a large range of metallicities. The model pre- 
dictions are from Francois et al. (2004) and are normalized to the predicted solar abundances. In the 
top right of each figure we indicate our predicted solar ratios, all in good agreement with the observed 
ones([a/Fe]=0.0). The black dots represent the data from Cayrel et al. (2004). For the other sources of 
data see Francois et al. (2004). 
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FIGURE 4. The same as in Figure 3 for other chemical species. The predictions are from Francois et 
al. (2004). The black dots are the data from Cayrel et al. (2004). For the other sources of data see Francois 
et al. (2004). 



stars of the Milky Way. Previous work on the subject had shown a large spread in the 
abundance ratios of these elements to iron, especially at low metallicities. This spread 
is confirmed although is less than before, and is at variance with the lack of spread 
observed in the other elements shown before. Apart from this problem, not yet solved, 
these diagrams can be very useful to put constraints on the nucleo synthetic origin of 
these elements. In particular, Cescutti et al. (2006) by adopting the two-infall model 
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FIGURE 5. The evolution of deuterium in the solar vicinity during the galactic lifetime. The assumed 
primordial D value is from WMAP (small bar at t=0), the large bar at t=0 represents the primordial D as 
measured from QSO absorption spectra. The bar at t=8.5 represents the solar system value from Gleiss 
& Gloecker (1988). The D abundance in the local ISM (at t=13 Gyr) is from Linsky (1998, small bar) 
and from FUSE and HST-GHRS data (large dotted bar, Vidal-Majar et al. 1998; Moos et al. 2002). The 
thick line represents the predictions of the two-infall model whereas the thin line the predictions of Tosi's 
(1988) model, as described in Romano et al. (2003). 



predicted the evolution of [Ba/Fe] and [Eu/Fe] versus [Fe/H], as shown in figures 8 and 
9. They can well fit the average trend but not the spread at very low metallicities since 
the model assumes instantaneous mixing. In order to fit the Ba evolution, they assumed 
that Ba is mainly produced as s-process element in low mass stars but that a fraction of 
Ba is also produced as an r-process element in stars with masses 10-30M Q . Europium is 
only an r-process element produced also in the range 10-30M Q . 



CONCLUSIONS 

Our main conclusion on the origin of the elements and the mechanism of formation of 
the Milky Way are: 

• The elements 12 C and 14 /V are mainly produced in low and intermediate mass stars 
(0.8 < M/Mq < 8). The amounts of primary and secondary N are still uncertain. 
Nitrogen should be primary in massive stars. 

• The a-elements originate in massive stars: the nucleosynthesis of O is rather well 
understood and there is agreement between different authors. 

• Magnesium is generally underproduced by standard nucleosynthesis models. Start- 
ing from the yields of WW95 we suggest that to fit the data the Mg production 
should be increased in stars with M < 20M Q and decreased in stars with M > 20M Q . 
Silicon should be only slightly increased in stars with M > 4OM (Francois et al. 
2004). 
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FIGURE 6. The evolution of Barium in the solar vicinity as predicted by the two-infall model (Cescutti 
et al. 2006). Data are from Francois et al. (2006) 
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FIGURE 7. The evolution of Europium in tthe solar vicinity (Cescutti et al. 2006). Data are from 
Francois et al. (2006) 



• Iron originates mostly in type la SNe, and the Fe production from massive stars is 
still uncertain. The metallicity dependent yields of Fe from WW95 overestimate Fe 
in stars with M < 30M Q . It is better to adopt the Fe yields for solar metallicity all 
over the whole metallicity range. The Fe yields of 199 do not show this problem 
and they are computed only for solar metallicity (Francois et al. 2004). 

• The Fe-peak elements: Cr, Mn, Ni, Co. Always starting from WW95 we suggest 
that the production of Cr and Mn should be increased in stars with 10-20 M Q . 
Cobalt production should be higher in type la SNe and lower in 10-2OM Q , Ni 
production should be lower in type la SNe. The type la SN yields adopted for 
comparison are those of 199 (Francois et al. 2004). 
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• In order to reproduce the evolution of Zinc relative to Iron and in particular its 
overabundance at very low metallicity, new yields from massive stars are required. 
They should produce more Zn than the standard yields. A good approximation is to 
use the Zn yields from WW95 for solar metallicity. Moreover, Zn production from 
type la SNe should also be increased (Francois et al. 2004). 

• Barium is mainly an s-process elements produced in the 1-3M© range, but it must 
have also an r-component from 10-30M© stars. Europium is only an r-process 
element formed in the stellar mass range 10-30M Q (Cescutti et al. 2006). 

• Deuterium is only destroyed in stars. A very good agreement with observations is 
obtained when the primordial abundance (D/H) p = 2.5 • 10~ 5 derived from WMAP 
results is adopted, and the astration factor during the galactic lifetime is < 1.5 
(Romano et al., 2003). 

• In order to fit the G-dwarf metallicity distribution in the solar vicinity, a timescale 
for the thin disk formation at the solar galactocentric distance of 7-8 Gyr is required. 
To reproduce the observed abundance gradients along the galactic disk an inside- 
out formation mechanism with the timescale for disk formation increasing with 
galactocentric distance should be assumed (Matteucci & Francois, 1989; Chiappini 
etal. 2001). 
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